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ABSTRACT 


Conformal mapping was used to obtain the free streamline pattern and surface pres- 
sure distribution for downward jet flow issuing from two parallel slot nozzles. After 
striking the ground, a portion of the flow from each jet moves along the ground toward 
the centerline between the jets. When these two portions collide, an upflow results 
which strikes a plate above the ground simulating the undersurface of an airplane fuse- 
lage. Typical flow patterns are shown to illustrate the effect on the free streamlines of 
nozzle height, nozzle spacing, plate height, plate width, and plate angle. 
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SUMMARY 

Conformal mapping was used to obtain the free streamline pattern and surface pres- 
sure distribution for downward jet flow issuing from two parallel slot nozzles. After 
striking the ground, a portion of the flow from each jet moves along the ground toward 
the centerline between the jets. When these two portions collide, an upflow results which 
strikes a plate above the ground simulating the undersurface of an airplane fuselage. 
Typical flow patterns are shown to illustrate the effect on the free streamlines of nozzle 
height, nozzle spacing, plate height, plate width, and plate angle. 


INTRODUCTION 

Under a vertical -takeoff (VTOL) airplane, as illustrated by figure 1, there is a com- 
plicated flow pattern. The downward directed jets in figure 1(b), react against the 
ground, and a portion of the flow is turned outward and flows along the ground. The 
remaining flow moves inward, and under the fuselage the opposing streams collide and 
move upward. Portions of this upward flow will move around the fuselage or will recir- 
culate under the wings. The upward flow in some instances provides a useful lifting 
force. However, it is generally undesirable because of recirculation of hot exhaust 
gases around portions of the airplane. This recirculating flow may be partially ingested 
into the engines, along with dust and debris from the ground that has been entrained 
(refs. 1 and 2). 

The complicated interaction with the ground and the airplane, and the presence in 
many instances of turbulent mixing, makes the analytical prediction of the three- 
dimensional flow flow field under a vertical -takeoff airplane extremely difficult. To 
obtain analytical solutions a simplified incompressible, isothermal flow model is employed 
herein that retains some of the major features of the actual flow field. The analysis is 





Figure 1. - Flow configuration for fan-pod VTOL aircraft. 


an extension of that in reference 3. As shown in figure 1, a fan -pod type of configuration 
is being considered. A two-dimensional approximation will be made to study the flow in 
the region between the engines. The flow is assumed to be inviscid; this could be approxi- 
mately true when the nozzle exit planes are within a few nozzle widths of the ground and 
within a few widths apart. For these conditions there will be only a small entrainment of 
surrounding fluid into the jet region prior to the flow being turned by the ground and 
turned under the fuselage. The two-dimensional model is shown in figure 2. A plate, 
which can be at an angle to the ground, has been used to simulate each half of the under- 
surface of the fuselage. 

By using a two-dimensional inviscid model, the determination of the flow field 
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Figure 2. - Two-dimensional model of jet deflection by ground and fuselage. 


becomes a free streamline problem in constant pressure surroundings. A solution can 
be obtained by using the Helmholtz -Kirchhoff method (ref. 4). This conformal mapping 
procedure is used to find a functional relation between the complex conjugate velocity and 
the complex potential of the flow. Then the flow configuration is obtained by integrating 
the complex potential multiplied by the reciprocal of the complex conjugate velocity. 

The flow pattern depends on the flow condition leaving the nozzle. The condition 
assumed herein is that the flow leaves perpendicular to the nozzle exit plane. This is a 
reasonable condition for nozzles of small width or when turning vanes are used to guide 
the flow. The interaction of the flow with the ground produces a nonuniform downward 
velocity leaving the nozzle; the velocity in the central region of the nozzle exit plane is 
less than that along the nozzle sides. 

The final analytical results were evaluated for several combinations of the param- 
eters, such as plate height, plate width, and spacing between nozzles. Typical free 
streamline flow patterns are given along with the pressure coefficient along the plate, and 
the velocity distribution across the nozzle exit plane. In addition to revealing the nature 
of the flow, the inviscid solution provides the upper limit of the pressures achieved on the 
ground and fuselage. The inviscid solution is also the zeroth order configuration which 
is needed to compute viscous entrainment by the flow. 


SYMBOLS 


A a constant 

A q a constant 

Cp pressure coefficient 
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quantity defined in eq. (A8) 

flow region in cp- plane 

flow region in W-plane 

complete elliptic integral of first kind 

defined by K'(k) = K(k') 

modulus of elliptic integral 

defined as |/l - k^ 
quantity defined in eq. (18) 
pressure 

ambient pressure outside jets 

functions defined in eqs. (B4) and (B5) of ref. 3 

functions defined in eqs. (B2) and (B3) of ref. 3 

complex variable in T-plane, £ + i-q 

complex variable in t-plane 

velocity in x-direction 

velocity 

velocity at nozzle exit 
average velocity at nozzle exit 
velocity along free streamlines 
velocity in y-direction 
complex potential, $ + iip 

coordinates of center of nozzle exit plane nondimensionalized by A 

coordinates in physical plane 

complex variable, x + iy 

tilt angle of plate, fig. 3 

flow region in T-plane 

angle of deflected jet, fig. 3 

width of nozzle 

widths of streams flowing to left and to right 


£ complex conjugate velocity, u - iv 

p density 

£, r] coordinates of T-plane 
0 quantity defined in eq. (A7) 

$ velocity potential 

x[/ stream function 

O function defined in appendix A of ref. 3 

Subscripts: 

A,G,H refer to stagnation points, fig. 3 
f fuselage 

g ground 

L flow to left 

R flow to right 


ANALYSIS 

The analysis presented herein is a generalization of, and quite similar to, that given 
in reference 3. For this reason only those parts of the analysis differing from reference 3 
will be discussed in detail. The configuration of the flow, as shown in figure 2, is sym- 
metric about a vertical plane between the nozzles, hence only half the flow field need be 
considered. The boundaries of the jet in the physical plane (z-plane) are shown in figure 3. 
Since the region outside the jet is at constant pressure, Bernoulli's equation shows that 
the velocity has constant magnitude along the free streamlines 6 e*, and The 

direction of the velocity must be along the boundaries on Xb', XoT, and fjHEl corresponding 
respectively to the bottom of the fuselage, the vertical line of symmetry, and the ground. 
The point H is a stagnation point at which the direction of the velocity along (jHET must 
reverse. The dashed line is the dividing streamline of the flows going to the right and 
left. The points G and A are also stagnation points. Along the line the velocity is 
vertically downward. 

As in reference 3, let £ be the complex conjugate velocity u - iv. From the spec- 
ified conditions on the velocity at the boundaries of the flow field, it can be deduced that 
the flow field in figure 3 maps into the interior of the region I of the hodograph plane 
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Figure 3. - Flow boundaries in physical z-plane (z * x + iy). 



Figure 4. - Flow region in hodograph £ -plane 
(£ - u - iv). 


shown in figure 4. Let W = $ + be the complex potential for the flow, where $ is 
the velocity potential and \jy is the stream function. Considerations analogous to those 
in reference 3 show that the flow field maps into the interior of the region J of the complex 
potential plane in figure 5. 

The solution in the physical plane is obtained from the integration 

z = J' — dW + constant (1) 

To carry out this integral, £ and W must be expressed in terms of the same vari- 
able of integration. This is done by appropriately mapping the rectangular region r of 
an intermediate T-plane (described in ref. 3 and depicted in fig. 7) into the regions I and 
J of the £ and W planes, respectively. The corresponding positions of the boundaries 
of r, I, and J are indicated by the lettering scheme in figures 4, 5, and 7. 


G ABF 

F 

Figure 5. - Complex potential plane (W = $ + i^). 
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Mapping Function Between £ and T-Planes 

The function £ which properly maps the rectangle r into the region I of the 
£-plane can be constructed by using the function Cl defined and studied in appendix A of 
reference 3. By a procedure completely analogous to that in reference 3, it can be 
demonstrated that the proper mapping function £ is given by 

£(T) = V o e i(7r/2 W; * H ,k)V«(T; %,k) [h(T; ^k)]^ 77 ; Ter (2) 

provided that the relation between £jj, and £ A is chosen to satisfy the following 
requirement. Along the plate the argument of £(T) must equal the direction along 
the plate, so that from figure 4 


arg £(T) = — + j3; T e && 

2 

It follows by the use of equation (A24) of appendix A of reference 3 that the argument 
of equation (2) is 


arg 


£(T) = i + i(^ + x \ + il^. + 1\ + £(i^ + A. 

2 2 \K ) 4 \K / 2 \K ) 


T e AB 


Equating the previous two expressions for arg £(T) gives the condition relating 
l G , and g A : 




( 3 ) 


Mapping Function Between W and T Planes 

The function W which properly maps the rectangle r of the T-plane (fig. 7) into 
the region J of the W-plane (fig. 5) can be constructed by use of the intermediate t-plane 
shown in figure 6. An application of the Schwarz -Christoff el transformation shows that 
the mapping that transforms the upper half t-plane onto the region J of the W-plane in 
the manner indicated in figures 5 and 6 is defined by 
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Figure 6. - Intermediate t-plane. 


dW 

dt 


t - t 


= iA_ 


H 


(l + kt)(t_ - 1) Yl - t 


Vi 


9mt > 0 


(4) 


Now it is shown in reference 3 that the rectangle r of the T-plane is mapped onto 
the upper half t-plane (fig. 6) by 

t = sn(T, k); Ter (5) 

This can be combined with equation (4) to eliminate t and find the desired mapping 
from the T-plane to the W-plane. The relation given in reference 3 is used 



1 

ksn(£ H ,k) 


along with the relation 


cn(i77 F ,k) 

t F = sn(T F , k) = sn(K + i r^, k) = 

dn(iry F , k) 


1 cn(?7 F , k' ) 

cn(? 7 p, k') dn(77 F ,k') 


1 

dn(?? F , k T ) 


Upon carrying out this procedure the required function that maps T onto J is found 
to be 

j ... (ksnf^snT - 1) dnT 

^W = iA % J.. Ter (6) 

dT (1 +ksnT) 1 - dn(? 7 F , k')snTj 


where A is a new constant equal to 

dn^k’) 
° ksn^ H 
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By using partial fractions equation (6) can be rewritten as 


dW iA J ^ sn ljj _ dn(?7p^k') k(l + sn£jj) 

d(snT) k + dn(r ?F ,k’) j / ~ j ~ 

T 1 - sn T [l - ^(Tjpjk'JsnTj Vl - sn^T (1 +ksnT) 


Integrating this expression and neglecting an unimportant integration constant give 



-k + iio 


K + iK' 



E 

H 

G 

A 

B 

FC 




D 

r 



C 



-K + iO 




K + iO 


Figure 7. - Mapping of jet flow into region r of T-plane (T = £, + \ri). 

The flow through the nozzle is equal to the difference between the stream function at 
points D and C. Noting that W = $ + ip and that $(C) = <3?(D), the flow must equal 
i//(D) - xjy( C) = |w(D) - W(C)]/i. When the coordinates of D and C in figure 7 are used, 
the flow through the nozzlebecomes: [w(-K + iO) - W(K + iO^/i. Hence, if Vj is the 

average velocity at the nozzle exit, 


V A = W (~K + iO) - W(K + iO) 
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Using the fact that sn(±K) = ±1 and choosing the appropriate branches for the inverse 
sine, we find from equation (7) 


W(K) = 


k’sn(rj F 


i[dn(rj F ,k’) -ksn^Jl - - k(l + sn^snfap, k’) ll 

,k')[k + dn(T 7 F ,k')] l L F ^2 H * 2j 


W(-K) = , <fdn(77 F ,k') -ksn£ H ~| — - k(l + sn£ H )sn(? 7 F ,k’)(-- 

k'sn(? 7 F ,k')[k + dn(ry F ,k')] 1~ * ^2 \ 2, 


Hence, 


V T A = p- -fdn(77 F ,k') 

1 k’snC^jkOfk + dn(r7 F ,k')J L 


- k sn^ H 


+ k(l + sn 



which gives the constant A as 


V T A 


A = — ? — k'sn(77 F ,k') <1 + (1 + sn£ H )k 


n ( 7 7p> k' ) ' 
dn(?7 F ,k') + k 



( 8 ) 


An alternate form for A will be given by equation (15). 

As shown in figure 3, the interaction of the jet with the ground causes it to divide. 
The asymptotic width of the flow region to the right is designated by 6 R , and that to the 
left by 6 l . The flow to the left must equal the jump in the imaginary part of W at the 
point F. In order to evaluate this jump, notice that F is at T = K + i?^ along the 
boundary T = K + irj in figure 7 and that 

sn(K + i? 7 ) = - 

dn(? 7 , k T ) 

Then, along this boundary, the inverse sine terms in W (eq. (7)) become 
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sin - * 


dn^k') _ 


dn ( 77 , k’) | 


d n{rijp, k') 
dn(i7,k») 


- 1 


= -i lnli 


jl - dn(?7 F , k') dn( 7 j,k’) + ^/[l - dn 2 (77, k' )] [l - dn 2 (?7 F , k’)j] 
dn(T7,k') - dn(?7 F ,k') 


= - - i lnff^)] + i lnjdn(?7, k') - dn^.k’)] 


where 


f 1 (7 7 ) - 1 - dn(T7 F ,k') dn(77,k') + ^[l - dn 2 (T7,k’)] [l - dn 2 (^ F ,k')] 


sin -1 


k + 1 


dn(77,k') 

-if 

1 k 

_ dn(T 7 ,k’)_ 


L k dn(77,k’) + f 
dn( 77 , k’) + k 


dn(7j,k') + k 


7T 

2 



where 


f (r? ) = k dn( 77 ,k’) + k’ 2 sn( 77 ,k’) + 1 
2 k + dn( 7 j, k’) 

Since 0 < dn(77, k’) < 1 and sn(7j, k’) is positive for 0 < 77 < K’ it is found that 
f 1 (7 ? ) > 0 and f 2 ( 77 ) > 0 for 0 < 77 < K 1 and, consequently, the logarithms of f^ and fg 
are real. 

Hence, choosing the appropriate branch for the lnjdn(77,k') - dn(7T F ,k')J term and 
substituting in equation ( 7 ) gives 
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II II 1 1 I II I IIH1II ■ 1 1 I III HIM II llll ■■ HUM II 


9m W(K + irj) = T ji fdn( 7 ? F , k’) - k ?u|J] 

k's n (?7 F , k')[k + dn^k')] |2 L HJ 


- | (1 + sn^sn^k’)!; v > Vp (9) 


JW(K + irj) = ^ {- i [dnjrfc , k’)-k sn|J| 

k'sn(77 F ,k’)[k + dn(77p,k-)] \ 2 L ^ HJ 


- — (1 + sn ^jj) sn ( ? 7pj k')j; r] < ? 7 f (10) 


Since the velocity at the point F must be V it follows from equations (9) and (10) that, 
if 6^ is the asymptotic perpendicular jet width (fig. 3) and e > 0, 

6 l V q = lim 9m[ w[k + i(7? F + e)| - w[k + ity, - e )j} 


Air 


k'sniriy, k’)[k + dn^k’)] 


7n - ksn^J 


An 


k'sn(77 F ,k') 


1 - 


k(l + sn£ R ) 


k + dn(?7p, k') 


Equation (8) is used to eliminate A giving 


6 L V o = V I A 


k + dn(r 7 F ,k’) - k(l + sn| R ) 


k + dn(? 7 F , k') - k(l + sn£ R )[l - sn( 7 j F ,k’)] 


(ID 


Continuity requirements dictate that 


5 p V n + 6 t V = AV t 
R O L O I 


Hence, using equation (11), the 6^ can be eliminated to obtain 


5 R V I 


k(l + sn ^pj)sn k T ) 


A V Q k + dn^k') - k(l + sn£ R )[l - sn( 7 j F ,k')] 


( 12 ) 
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From equations (12) and (11) the ratio of asymptotic widths is 


6 r _ sn(r 7F ,k’) 
S L k + dn(T ?F ,k') 
k(l + sn^ H ) 


By rearrangement 


k(l + sn£ H ) 

k + dn(77p,k') 


1 


6 I 

— sn(77 F ,k’) + 1 

6 R 


(13) 


and 


— dn(? 7 F ,k') -ksn(rj F ,k') 

6 t 

sn % = ^ 

— k +ksn(77 F ,k') 

6 L 


Using equation (14) in equation (8) gives the constant A as 


Vi A 

A = — - — k’ 


sn(r7 F ,k’) + — 
6 L 


1 + 


R 


(14) 


(15) 


Integration to Yield Coordinates in Physical Plane 

When the complex conjugate velocity £ and the complex potential W as functions of 
the parametric variable T are known, the physical variable z can be found as a function 
of T by using equation (1). The integration is carried out in the T-plane instead of in 
the W-plane, so that equation (1) becomes 
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-r -i 

Jij r +iK* 


z = , ^dT 

£(T) dT 

5 G ' r ' 


(16) 


The origin of the coordinate system has been chosen at the point G. Substituting equa- 
tions (2), (6), and (15) in equation (16) gives 


z 

A 


J_\ 

v o h 


sn(77 F ,k>) + — 
6 L 


1 + 




X 


f 


(ksn|jjSnT - 1)(1 -ksnT) dT 


f dnT[l - dn(? ?F ,k’)snT]s2(T; | H ,k)^(T; £ G ,k)[f2(T; 


(17) 


Average Velocity at Nozzle Exit 

To find an expression for Vj/V Q , note that the nozzle width is given by 

A = z(-K + iO) - z(K + iO) 

Then it follows from equation (17) 


V I 

sn(?] F ,k’) + 

V o 

6 r 

1 + 3 

L 6 l 


U R 


■«-f f -r 

r a dns E 


,-l 


(1 -ksn£ H sn£)(l -ksn|)d| 


(18) 


- dn(?7p, k’)sn^jsi(^; ? G ,k) jft(£;? A ,k)j 


P/ n I 


The quantity M defined in this equation will be used subsequently. 

Pressure Coefficients 

The pressure coefficients C and C f along the ground and fuselage are defined 

Pj & P? * 

by 
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P(x, 0) - p o 



C 


P,f 


p(Xf,y f ) - p q 



where Xf,y f are the coordinates of the points along the fuselage. Using Bernoulli’s 
equation gives 


C P,g 1 


u(^K’) 


n2 


V„ 


- K - £ — £f 


C . = 1 - 
p> * 


u 2 (^K’) + v 2 (g,K') 




(19) 


( 20 ) 


Summary of Analytical Relations 


For convenience, the results of this section are now collected in one place. 


6 r 

— dn(77 F , k ’) -ksn(? 7 F ,k') 

§T 

sn ^H g — (14) 

— k + ksn( 7 v.,k') 

5 t r 


I + £ 4K = M 

4 2/ n n u 77 A 


(3) 


= iO(T; | R ,k) yj2(T; i^k)[n(T; l A ,k^ /,r (2) 

* /-» 
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1 


W(T) _ i 


Vj A it [dn( 77 F ,k’)+k] + (1 + sn£ H )k [sn^jk’) - l] 

Ijdn^,^) -ksn^J sin * 


^(^Tpjk') - snT 


dn (’?F» k ')snT - 1 

- k(l + sn^ H )sn(77p,k') sin" 1 | " k + snT ~ |l 
w * U+ksnTj 


z(T) -( l\ k ’ 

A ~Voh 


sn(77 F , k ’) + — 

5 L 


1 + 


R 


(ksn^snT - 1)(1 -ksnT)dT 


? dnTjl - <111(7^, k')snT|o(T; | R ,k) fa(T; £ G ,k)[fi(T; * A ,kj[ 


P/ir 


6 R V I 


k(l + sn| H )sn(77 F ,k’) 


A V Q k + dn(j7 F ,k') - k(l + sn^fl - sn( 7 j F ,k’)] 



r 6 r " 

sn 0? F , k ’) + — 

V I 

6 l 

V o 

6 r 

1 + — 


L 6 l J 


=»= r T 

r L *‘E 


(1 - k sn^ H sn^)(l - k sn£) d£ 


- dn(7? F ,k')sn^(4; { H , k)j 0($; ? G , k) |n(fc ( A , k^ /n I 


C p,g = 1 - 


uU>K’) 


C f = l 
P, f 


u 2 (^> K*) + v 2 (|, K*) 


‘ K * * * «G 
; * € * K 


Additional forms of equations (2) and (17) that are useful for computer evaluation are 
given in the appendix. 


(7,8) 

(17) 

( 12 ) 

(18) 

(19) 

(20) 


16 



COMPUTATIONAL PROCEDURE 


If all coordinates are nondimensionalized by the nozzle width A, there are five 
independent parameters governing the flow configuration: plate (fuselage) height above the 
ground, plate width, plate angle, nozzle height, and spacing between nozzles. The plate 
angle /3 appears explicitly in the analytical expressions and can be directly specified as 
an input variable in the computer program used to evaluate numerical results. The other 
four physical quantities must be found by computation in terms of four convenient input 
quantities all having values between 0 and 1: k, 6 L /S R , 77-p/K', and K. From a 
chosen k, the k', K, and K' can be found from 



K'(k) = K(k') 

Then can be found from equation (14), and all quantities which are necessary to 
determine from equation (3) will then be known. The quantities M and Vj/V Q are 
then found by carrying out the integration in equation (18), using the S2 functions eval- 
uated by the method described in the appendixes of reference 3. The dimensionless width 
of the stream flowing to the right 6^/ A is evaluated from equation (12), and 6^/ A is 
computed from (6 j,/a)(5j_ / 6^), where 5-^/6^ is one of the specified input quantities. 

With all these quantities evaluated, the height of the plate is found from equation (A6), 
and its width by integrating equation (A5) to an upper limit of K. These dimensions, 
along with the specified angle / 3 , fix the position of point B. The streamline 131* can 
then be plotted by use of equation (A12). The distance in equation (A10) is then evaluated 
to fix the horizontal position of the nozzle. The vertical separation of points D and E 
along the right free streamline is found from the second of equations (A9) by integrating 
to an upper limit of K* . This distance, along with 6p/A, is used to determine the noz- 
zle height. The free streamlines originating at the nozzle are then computed from equa- 
tions (A9) and (All) and are drawn starting, respectively, from points D and C. 

The pressure coefficients along the ground and fuselage are found from equations (19) 
and (20) using the velocities evaluated from the second equations of (A4) and (A5). 
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RESULTS AND DISCUSSION 


As revealed by the analysis, there are several independent parameters (all made 
dimensionless by dividing by the nozzle width) governing the flow: the nozzle height, 
spacing between nozzles, plate height, plate width, and plate angle. It was not feasible 
to systematically compute flow patterns for the wide variety of possible combinations of 
these parameters. To limit the number of computed flow patterns somewhat, the plate 
simulating the underside of the fuselage was fixed in a horizontal orientation with the 
exception of the results shown in figure 12 where a comparison is made with a plate tilted 
at 45°. The plate was generally positioned at the same height as the nozzle, and only two 
nozzle heights were considered (Y^ = 1 and 2). 

Consider figure 8(a) as a typical set of results. On each set two cases are given, one 
in solid lines and the other dashed. Each case shows the velocity distribution across the 
nozzle exit plane, the free streamlines bounding the moving fluid, the pressure coefficient 
along the ground, and the pressure coefficient along the plate (fuselage). 

Figures 8(a) and (b) illustrate the effect of moving the nozzle and plate upward to- 
gether, as when the airplane is taking off. Part (b) is for a longer plate than part (a), but 
both portions of the figure exhibit the same general features. In the upper position the 
exit velocity from the nozzle is more uniform; hence, the flow is increased as the nozzle is 
raised since there is less reaction from the ground. The increase in flow results in the 
integral of the ground pressure coefficient being a little larger for the upper nozzle posi- 
tion, thus providing a small increase in lift. The pressure coefficient acting on the under- 
side of the plate is decreased a little as the plate and nozzle are raised, thereby reducing 
the lift on the fuselage. In figure 8(b) the stream deflection by the plate is so large for 
the upper nozzle position (dashed) that the stream passes back into the region of the nozzle 
exit. In an actual flow the two streams would collide. The analysis used herein does not 
allow for such interference and permits the two streams to pass independently through 
each other (on different sheets of the complex plane). Actually, there would be interfer- 
ence and probably a recirculating flow region for this configuration. 

Figure 9 demonstrates the effect of changing the horizontal position of the nozzle. 

Parts (a) and (b) are for a low nozzle (Y N = 1), while parts (c) and (d) are for a high noz- 
zle (Y n = 2). The dominant effect is that the flow going to the left and impinging under 
the fuselage becomes quite small as the nozzles approach the plate. Thus having the fan 
pod adjacent to the fuselage rather than as shown in figure 1(a) will reduce the hot gas 
circulation around the fuselage. Recall that this is a two-dimensional analysis. If the 
nozzles were round rather than in a slot or pod configuration, there would be an additional 
effect of nozzle spacing resulting from the radial spreading of the flow as it moves out- 
ward along the ground. This effect along with entrainment would tend to decrease the flow 
under the fuselage as the spacing between nozzles is increased. 
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Dimensionless coordinate, x/A 


(b) Low nozzle (Y N = 1) ; plate width -0.8. 
Figure 9. -Continued. 
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In each of figures 10(a) and (b), the plate width is varied while the nozzle is at a 
fixed position. The (a) and (b) parts are for nozzle heights Y N = 1 and 2, respectively. 
For the wider plate the deflection of the left -moving stream is increased and there is a 
large stagnation region below the plate, especially when the nozzle is in the lower posi- 
tion. As the plate width is increased, the space between the nozzle and plate is dimin- 
ished, thereby reducing the flow to the left. Raising the plate while keeping the nozzle 
fixed also increases the jet deflection, as shown in figure 11. 

Figure 12 shows the effect of tilting the plate while keeping fixed the horizontal pro- 
jection of the plate and the nozzle position. The flow pattern is changed only a small 
amount, the stream flowing to the left being turned more by the tilted plate. The exact 
shape of the underside of the fuselage is thus of minor importance. 

The previous figures have given the reader a quantitative appreciation of the flow 
trends resulting from varying the five independent parameters governing the flow config- 
uration. The flow patterns can also be used as the zeroth order solution for evaluating 
entrainment into the flow regions. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 24, 1969, 

129-01-07-07-22. 
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APPENDIX - WORKING FORMULAS FOR COMPUTING FLOW PATTERN 


In this appendix equations (2) and (17) will be rewritten in various forms which are 
convenient for computing the velocities on and the positions of the various boundaries of 
the flow. To accomplish this the quantities defined in appendix B of reference 3 will be 
used. The positions along the free streamlines are found by integrating equation (17) 
along the sides of the rectangle of figure 7 where T = i 77 ±K. By using the various for- 
mulas for the change of argument of elliptic functions, the following relation is obtained 
in which the elliptic functions depend on real arguments: 


/dw\ = ivI A k’ 2 

6 R 

sn(7j F ,k')+ — 
6 L 

|±k sn(£ H , k) - dn(7j,k')]cn(r7,k') 

\ dT /T=ir?±K 77 

r 

+ 

O) I Oj 

f 1 5d 

i 

|dn(? 7 ,k') ±kj|dn(? 7 , k') t dn^k')] 


Similarly, along the top of the rectangle, which corresponds to the solid boundaries and 
plane of symmetry, 



T=£+iK' 


Vj A k’k 

6 r ” 

sn(77 F ,k’)+ — 
6 L 

[sn(| H ,k) - sn(£, kj|cn(£, k) 

7 r 

6 r 

1 + — 

6 l J 

[l + sn(4,k)][ksn(|,k) - dn^,: 


(A2) 


Using the results and definitions of appendix B of reference 3 and substituting equa- 
tions (Al) and (A2) in equations (2) and (17) give the following working formulas (note 
that M was defined in equation (18)): 

Position and Corresponding Velocity Along Nozzle Exit 


x(-K, 0) - x(|, 0) _ X D ~ X 


M : 


/ 


vM) 


[l - k sn(£ H , k)sn(£, k)J [l -ksn(£, k)] d£ 


dn(|,k)[l - dn(i^,k')snU,kj|Rg\0 [kjf >(« jf 7 * 

= f^ 4 ) u) 


^ -K < ^ < K 


(A3) 
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Position and Corresponding Velocity Along Ground 


x(g,K') _ x ~ X G 


*\ 


Mkk' 


/ 

J ki 


[sn(£,k) - sn(? H ,k)][l - sn(|,kj] dj 


cn(«,k)[dn( 7 F ,k-) -ksn^kj^d) ^ 


> -K s 6 * «g (A4) 


= r|j ) (?> V R G 1) <*>l?i 1) <*il 


P/tt 


Position and Corresponding Velocity Along Fuselage 


Z(5 + iK ') - iy ( ^ A . K ') = - Mkk ' 

A V 


t 


[sn(*,k) - sn(4 H ,kjj[l - sn(i|,kj] d? 


cn(^,k)[dn(77 F ,k') - k sn(£, k)]R^(?) V‘ R G )( ^[' R A^i] 


l (A5) 


|CU jjKMl - 

v„ 


u 2 (?, K ’) + v 2 ( 4 , K ') 
V o . 


1/2 


-R^tt) V- R G )(?) [- R A )( ^ 


/3/7T 


Height of Point A on Fuselage Above Ground 


y(H a ,k') y A - y G -Mkk- 


*«A 


[sn(|,k) - sn(| H ,k)][l - sn(£,k)] d£ 


cn($, k)(dn(? 7 p,k') -kan^k)]^ 1 ^) ^ -R^^R^d)]' 


Ph 


(A6) 
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Coordinates of Free Streamlines 


Define and H ± as follows 


0 (i) (7?) = - + q£\v) + - Q ( g ] (v) + £ Q i = 3, 4 

2 2 n 


± _ [dn(??,k') ±ksn(^ H ,k)]cn(77,k’) 

^ [dn(? 7 , k') =f kjfdnfok’) idn^,^)] 


(A7) 


(A8) 


Then the coordinates along the the right free streamline extending between the 
points D and E are 


X(-K, „) - x(-K, 0) _ k^M /"> H+(lj) cos r e (S) ( 3 dl) 

A A 7T Jq L J 




> 0 < 77 < K' 


y(-K, „) - y(-K, 0) , . . k^M /”> H+(j)) sln foP), H dl) 

a a 7r y 0 L J 


(A9) 


x(K,K')-x(K,0) t X B - X C = M py J*' H -(„) cos [e< 4> (r,)] d, - sm [® (4> <7^.)] 


(A10) 


The last term is the horizontal displacement of the left branch of the jet at infinity. 
The coordinates along the free streamline between C and F are 


x(K, ij) - x(K, 0 ) _ x ~ X C _ k'*M /•” H -(,) cos |e (4) („)] d„ 

A A 7T Jq L J 

y(K, v) - y(K, 0 ) = y l y C = _ kl 2 M f V B - {v) sin [0^(77)] d7? 
A A it J n L J 


> 0 < 7? < 7 ?p 


(All) 


The coordinates along the free streamline between B and F are 
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I 


x<K t K') -x<K t i) . ^ = id M f H- (,) cos [e< 4 >(„)] d„ ' 
y<K,K-)-y(K, n ) m ] d M y K ' H-(,) sin [e< 4 ><„)] d„ 


7f < v ^ K ’ 


(A12) 


33 



REFERENCES 


1. Spooner, S. H. : The V/STOL Aircraft Environment. Paper No. 68-GT-40, ASME, 

Mar. 1968. 

2. Kemp, E. D. G. : Studies of Exhaust Gas Recirculation for VTOL Aircraft. Paper 

No. 67-439, AIAA, July 1967. 

3. Goldstein, Marvin E.; and Siegel, Robert: Two Dimensional Inviscid Jet Flow From 

Two Nozzles at an Angle to a Plane Surface. NASA TN D-5064, 1969. 

4. Birkhoff, Garrett; and Zarantonello, E. H. : Jets, Wakes, and Cavities. Academic 

Press, 1957. 


34 


NASA- Langley, 1969 1 E-4903 



